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ABSTRACT 

We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution 
iV-body simulation of the ACDM cosmology. We find that the cumulative mass distribution of specific 
angular momentum j in a halo of mass is well fit by a universal function, M{< j) = M^fxj /{jo + j)- 
This profile is defined by one shape parameter (/i or jo) in addition to the global spin parameter A. It 
follows a power-law M{< j) oc j over most of the mass, and flattens at large j, with the flattening more 
pronounced for small values of fj, (or large jo at a fixed A) . Compared to a uniform sphere in solid-body 
rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. 
High-A halos tend to have high fi values, corresponding to a narrower, more uniform j distribution. The 
spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each 
halo for ^ 80% of the halos. The more misaligned halos tend to have low-/i values. When averaged 
over spherical shells encompassing mass M, the halo j profiles are fit by j{M) oc M** with s = 1.3 ± 0.3. 
We investigate two ideas for the origin of this profile. The first is based on a revised version of linear 
tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on 
j transport in minor mergers. 

Finally, we briefly explore implications of the M{< j) profile on the formation of galactic disks assuming 
that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from 
an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep 
central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent 
of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed 
perfect j conservation. A possible solution is to associate the central excesses with bulge components 
and the outer regions with extended gaseous disks. 

Subject headings: cosmology — dark matter — galaxies; formation — galaxies: structure 



X 



1. INTRODUCTION 

The origin of the distribution of mass and angular mo- 
mentum in disk galaxies is still an open issue, despite its 
long history. 

The archetypical model of Eggen, Lynden-Bell & 
Sandage (1962) provides a useful framework, in which a 
disk galaxy (the Milky Way) forms by centrifugal support 
of a collapsing gas cloud. Mestel (1963) added the assump- 
tion that the specific angular momentum of each mass el- 
ement, j, is conserved during the collapse, and demon- 
strated that the final mass profile of the disk could then 
be related to the initial mass and angular momentum dis- 
tribution of the gas. The key is that the cumulative mass 
with specific angular momentum less than j, M(< j), is 
preserved during the collapse. Crampin & Hoyle (1964) 
(also Innanen 1966; Freedman 1970) then realized that 
the observed exponential density profiles of disk galaxies 
in circular motions are consistent with the angular mo- 
mentum distribution of a hypothetical uniform sphere in 
solid body rotation, 



M{< j) 
° Hubble Fellow 



1-(1-J7jmax)'/'' 



(1) 



where jmax = ^^R^ with u the fixed angular velocity and 
R the radius of the sphere, and Mtot the total mass of the 
sphere. 

Fall & Efstathiou (1980) re-examined this question in 
the context of the more modern view according to which 
disk galaxies form from contracting gas within extended 
dark-matter halos (White & Rees 1978). They added the 
assumption that the gas and dark matter were initially 
well-mixed such that the distribution of j in the disk is 
equal to that of the halo. Rather than assuming an an- 
gular momentum profile, they implicitly deduced it from 
the constraint that the final surface density profile of the 
disk is exponential. The implied disk scale radius, ra, is 
related to the total pre-coUapse angular momentum of the 
halo, roughly r^ oc A. Here, A is the dimensionless spin 
parameter (e.g. Peebles 1969) 

where J,E and M are the total angular momentum, en- 
ergy and mass of the system, and G is Newton's constant. 

Angular momentum is presumably acquired by the dark 
matter (and gas) through tidal interactions with neighbor- 
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ing objects (Peebles 1969). Barnes & Efstathiou (1987) 
(following Efstathiou & Jones 1979; Efstathiou & Barnes 
1983; Zel'dovich & Novikov 1983) utilized cosmological N- 
body simulations with a CDM power spectrum to make 
detailed predictions for halo angular momentum structure 
and statistics. In particular they found that the distri- 
bution of spin parameters is roughly log normal, with a 
median value of A ~ 0.05. They also confirmed that the 
angular momentum of pre-coUapse halos grows roughly lin- 
early with time, as predicted by linear tidal-torque theory 
(Doroshkevich 1970; White 1984). Several numerical and 
analtyical investigations have followed, which confirmed 
the nature of the A distribution, addressed the alignment 
of the spin vectors between neighboring galaxies, the (lack 
of) dependence of halo spin on environment, and the ef- 
fect of major mergers on halo spin (e.g., Frenk et al. 1988; 
Heavens & Peacock 1988; Zurek, Quinn & Salmon 1988; 
Warren et al. 1992; Catelan & Thcuns 1996; Cole & Lacey 
1996; Lcmson k Kauffmann 1997; Gardner 2000). 

Armed with this global information about the initial 
state of the dark matter and gas mixture, several authors 
have used the above disk formation scenario to predict ob- 
servable quantities. Blumcnthal et al. (1986) addressed 
the rotation curves of disk galaxies in the context of the 
CDM hierarchical formation scenario, assuming 10% of 
the mass in the gas component (Blumenthal et al. 1984). 
They pointed out that the non-dissipative dark halo should 
react to the dissipativc gas infall by considerable contrac- 
tion, and showed that this process can be approximated 
using an adiabatic invariant. This coupling between the 
halo and the disk "conspires" to produce a continuity be- 
tween the disk-dominated and halo-dominated regions of 
the rotation curves, as observed (Burstein & Rubin 1985). 

Flores et al. (1993) explored the dependence of the fi- 
nal rotation curves on a range of assumed parameters, in- 
cluding the measured spread in A. Dalcanton, Spergel, & 
Summers (1997) included a more realistic initial halo den- 
sity profile (Hernquist 1990) and explored the properties 
of both high surface brightness and low surface brightness 
galaxies. Mo, Mao, fc White (1998a; 1998b, 1999) per- 
formed a similar examination within the context of var- 
ious cosmologies, assuming that the dark halos followed 
the density profile advocated by Navarro, Frenk, & White 
(1996, 1997; NFW) and compared their predicted galaxy 
properties to both local and high-redshift data, van den 
Bosch (1998; 2000) and van den Bosch & Dalcanton (2000) 
extended the approach to include a bulge component and 
investigated the origin of the Hubble sequence of galaxy 
types and other aspects of galaxy formation. 

The main commonality in all of these investigations is 
that the variation in observable galaxy properties depends 
almost exclusively on the initial total spin parameter of the 
halo hosting the galaxy. However, all of the predictions 
rest on the implicit assumption that, if specific angular 
momentum is conserved, the initial angular momentum 
profiles of the halos produce disks that are exponential in 
form. ^ 

Within this simple framework for disk formation, the 
key missing ingredient is the actual distribution of angu- 

^ An exception is Dalcanton et al. (1997), who assumed that the initial 
disk density profiles that are similaj: to exponential (Ryden & Gunn 
Blumenthal, 1991). 



lar momentum in halos. Partially motivated by this, we 
set out to measure the j profiles directly using a statisti- 
cal sample of halos in a high resolution simulation, and to 
estimate the intrinsic scatter in these distributions. The 
long-term goal is to learn what the detailed form of halo 
angular momentum profiles can teach us about the origin 
of galaxy types, disk surface brightness profiles, and bulge 
properties. 

An additional strong motivation for this study comes 
from the results of hydrodynamic simulations of disk 
galaxy formation in a cosmological context. Unlike the 
simplified, monolithic- collapse framework discussed above, 
the hydrodynamic treatments model gas cooling within 
the hierarchical growth of dark halos, and typically find 
that the resulting disks are significantly smaller than real 
galactic disks (Navarro & Steinmetz 2000; Weil, Eke, & Ef- 
stathiou 1998; Navarro & Steinmetz 1997; Navarro, Frenk, 
& White 1995). The problem apparently arises because 
most of the mass is accreted through mergers with objects 
whose gas component has already cooled and clumped, 
causing a large fraction of their angular momentum to 
be transferred to the dark halo or transported outwards. 
Since angular-momentum transport is likely to always be 
from the inside out, the simple case in which j is conserved 
during the collapse, which we can study in detail with high 
accuracy, can serve as a useful limiting case. This will pro- 
vide a better basis for addressing the more complex phys- 
ical processes required for solving the angular-momentum 
problem seen in the hydro simulations. 

In §2 we describe our methods, including the simulation 
and halo finding, the measurement of the spin parameter, 
the associated errors, and the measurement of the angular- 
momentum profile. In §3 we discuss the measurcid (univer- 
sal) mass distribution of angular momentum in halos. In 
§4 we address the spatial distribution of angular momen- 
tum, its alignment, cylindrical symmetry, and profile in 
spherical shells. In §5 we use both tidal-torque theory and 
an estimate of angular momentum transfer due to minor 
mergers to explore the origin of the angular-momentum 
profile. In §6 we briefiy describe the implications for the 
disk surface density profile and the "angular momentum 
catastrophe". We discuss our results and conclude in §7. 
An Appendix is devoted to testing our results with a sim- 
ulation having 8 times the mass resolution, and investigat- 
ing M{< j) profiles of halos at high redshift. 

2. METHOD 

2.1. Halos in the ART simulation 

Large cosmological N-body simulations have reached the 
stage where detailed structural properties of many dark- 
matter halos can be resolved simultaneously. A method 
which allows the required force and mass resolution is 
the Adaptive Refinement Tree (ART) method (Kravtsov, 
Klypin, & Khokhlov 1997) which implements successive 
refinements of the spatial grid and time steps in high den- 
sity environments. The simulations based on the ART 
code provide a statistical sample of halos, covering a wide 
mass range, with the resolution appropriate for studying 
angular-momentum profiles. 

< j) profile was that of a uniform solid body rotator, which produces 
7) but only over about two disk scale lengths (Olivier, Primack, & 
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We have used the ART code to simulate the evolution 
of collisionless dark matter within the "standard" low- 
density flat ACDM model {fl,^ = 0.3, f^A = 0.7, h = 0.7, 
and (78 = 1.0 at z = 0). The simulation followed the trajec- 
tories of 256'^ particles within a periodic box of comoving 
size 60 ft.~^Mpc from redshift z = 40 to the present. A ba- 
sic 512'^ uniform grid was used, and six refinement levels 
were introduced in the regions of highest density, imply- 
ing a dynamic range of ~ 32,000. The formal resolution 
of the simulation is thus /res — 1.8/i~^kpc, and the mass 
per particle is nip = 1.1 x 10^ /i^^M©. 

Our halo finding algorithm (Bullock et al. 2000) is based 
on the Bound Density Maxima technique (Klypin & Holtz- 
man 1997) and has been specifically designed to identify 
halos (and subhalos) in such a high resolution simulation. 
In the current study of angular momentum, we limit our- 
selves to 'distinct' halos, which do not reside within a 
larger host halo, in the mass range 10^^ — 10^'' h~^MQ. 
In the Appendix, we use a higher resolution simulation to 
extend our mass range down to - 10" H-'^Mq. 

A brief description of the halo finder is as follows: Af- 
ter finding all the maxima in the smoothed density field 
of the simulation, we unify overlapping maxima, define a 
minimum number of particles per halo (say 50), and iter- 
atively find the center of mass of a sphere about each of 
the remaining maxima. We compute the spherical density 
profile about each center and identify the halo virial ra- 
dius i?v inside which the mean overdensity has dropped to 
a value Av, based on the top- hat spherical infall model. 
For the family of fiat cosmologies (fini + ilA = 1), the value 
of Av can be approximated by (Bryan & Norman 1998) 
Av ~ (IStt^ -H82a;- 39a;2)/(l + a;), where x = ^2™(z) - 1. 
In the ACDM model used in the current paper, Av varies 
from about 180 at z > 1 to Av ~ 340 at z = 0. 

We then fit each halo density profile with a universal 
functional form. We adopt the NEW profile. 




Fig. 1. — The distribution of halo spin parameter A' (histogram), 
compared to a log-normal distribution about Aq = 0.035 it 0.005 
with a width cr = 0.50 ± 0.03. 



2.2. Angular-momentum errors 

A first step towards measuring the distribution of spe- 
cific angular momentum within each of the halos is mea- 
suring the global spin parameter of each halo. The angular 
momentum of a halo of TV particles is defined by 



N 
i=l 



r: X V; 



(4) 



PNFwir) 



(r/rs) (1 + r/rs) 



(3) 



where rj and Vi are the position and velocity of the ith par- 
ticle with respect to the halo center of mass. In principle, 
given the spin parameter A, the value of the global specific 
angular momentum, J/M, can be determined by using an 
assumed energy content for the halo in Eq. (2). In prac- 
tice, however, this is not a straightforward procedure. For 
example, the energy of a halo in a crowded region is some- 
what ambiguous because it depends on the environment. 
A related difficulty arises when only a sub-volume of the 
virial sphere is concerned, e.g., within the cooling radius. 
We therefore define an alternative and more practical spin 
parameter by 

A' EE , (5) 

V2MVR 

given the angular momentum J inside a sphere of radius 
R containing mass M, and where V is the halo circular 
velocity at radius R, = GM/R. This spin parameter 
reduces to the standard A when measured at the virial ra- 
dius of a truncated singular isothermal halo. The value 
of A' turns out to be robust to the choice of an outer ra- 
dius; we use below the virial radius i?v, but we find a very 
similar distribution of A' values when using i?v/2 instead. 
This allows modelers to freely scale our results using any 
desired outer radius appropriate for the problem at hand. 

Removal of unbound particles is a new feature allowed by the high resolution of our simulation; it could not have been properly implemented 
in earlier studies. 



with the two free parameters and ps- An equivalent 
pair of parameters is, for example, the virial mass and 
the concentration parameter Cv = R^/rs- Using this fit, 
we iteratively remove unbound particles from each mod- 
eled halo, and unify every two halos that overlap in their 
rg and are gravitationally bound. ^ The modeling of the 
halos with a given functional form allows us to assign to 
them characteristics such as a virial mass and radius, and 
to estimate sensible errors for these quantities. 

The halo finding is 100% complete for halos of more 
than ~ 150 particles, M > 1.5 x 10" H'^Mq (see Sigad 
et al. 2000). For our purpose here, we limit the sample to 
halos of more than ^ 1000 particles. 
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Fig. 2. — Halo spin parameter A' versus concentration Cv. The 
symbol size is inversely proportional to the error on Cv. Shown are 
the linear regression line of A' on Cv (solid), and the relation ex- 
pected assuming that the standard spin parameter A is independent 
of Cv (dot-dashed). 



The spin parameters A' and A are in fact very simi- 
lar for typical NFW halos. They are related by A' — 
\\EsisIEnfw\^''^ ^ A/(cv)"^/^, where Esis and Enfw 
are the energies of virialized halos with isothermal and 
NFW density profiles respectively, and /(cv) ~ [2/3 -I- 
(cv/21.5)°-^] (Mo et al. 1998) ^. The function /(cv) is 
about unity for typical concentrations, Cv ~ 10. 

The distribution of A' over the ~ 500 halos in our sample 
is shown in Fig. 1. It is well fit by a log-normal distribu- 
tion. 



P(A') = 



1 



AV27rCT 



exp 



ln^(A7A[,) 
2a2 



(6) 
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with best fit values AJ, = 0.035 ± 0.005 and cr = 0.5 ± 0.3. 
Not surprisingly, the distribution of A' is very similar to 
the known distribution of A (Barnes & Efstathiou 1987). 
The distribution of A values for our halos has best fit values 
Ao = 0.042 ± 0.006 and cr = 0.5 ± 0.35. 

To highlight a small difference between A' and A, Fig. 2 
displays the values of A' versus Cv for our sample of simu- 
lated halos, and shows a weak correlation. Given that A is 
known to be uncorrelated with Cv (NFW), this correlation 
arises from /(cv) in the above relation between the two 
spin parameters (as demonstrated in the figure) 

The exact expression is /(cv) = 0.5cv[(l -I- Cv)^ — 1 — 2(1 -|- Cv) ln(l + Cv 
* In agreement with the NFW result, we also find no correlation between 



100 



Fig. 3. — Estimate for the error in J. The average (solid line) 
and standard deviation (dashed line) of J/Ji as a function of the 
number of particles in the random subsamples. 

The error on the measurement of J in the simulations is 
of primary importance for the interpretation of the results 
and, here, for the specific task of fitting a functional form 
to the angular-momentum profile. Obtaining a precise 
measurement of J in halos is difficult because the signal- 
to-noise is relatively small due to the fact that the local 
coherent velocity of the halo contributing to the angular 
momentum is small compared the velocity dispersion. The 
variable A' itself roughly characterizes the signal to noise 
(per particle) at any radius i?, Eq. (5), which is typically 
at the level of a few percent. 

In order to obtain a crude estimate of the error, which 
is probably an upper limit, we singled out several large 
halos with Ni simulated particles (> 2 x 10^), and mea- 
sured their angular momentum using a series of sparsely 
sampled realizations of iV particles down to 1% of Ni, with 
100 random realizations for each value of N. The average 
and standard deviation of the measured value J relative 
to the original value obtained with all the particles, Ji, as 
a function of N, are shown in Fig. 3 for three halos with 
different values of A'. We find that for a subsample of N 
particles the standard deviation aj scales like 

'^J /. i^r . i/(25A'2Ar) ^ J^. (7) 



J 



AVA^ 

This scaling relation is shown as error bars in Fig. 3, and 
it fits rather well the measured scatter. The first term in 
Eq. (7) corresponds to the general Poisson sampling error 
and the second term reflects the noise due to velocity dis- 
persion. The second term typically dominates because A' 
is small, and thus the spin parameter itself directly influ- 
ences our ability to measure J accurately, with the number 
of particles needed for a fixed accuracy scaling like A'~^. 
For A' = 0.03, even ^ 50% accuracy requires N 150. 
The above estimate of the error probably tends to over- 

)]/[cv-(l + Cv)ln(l + Cv)]2. 
Cv and the traditional spin parameter A. 
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estimate the true error because it includes small-scale noise 
that is not present when a smaller number of particles is 
used in the actual dynamical simulation of the halo. On 
the other hand, our procedure above does not mimic the 
error made in following the true dynamics of the smooth 
halo when a smaller number of particles is used — an ef- 
fect which tends to cause an underestimate of the error. 
An accurate evaluation of the error in J requires a direct 
comparison of the same halo as simulated with different 
resolutions, and is clearly sensitive to the A^-body code 
used. In our analysis here we limit ourself to the crude es- 
timate, Eq. (7). In the Appendix we use a simulation with 
8 times the mass resolution in order to check for systematic 
differences in the derived M(< j) profiles, and find none. 
This analysis gives us confidence that mass resolution does 
not affect our results in any severe way. 

2.3. Measuring the angular-momentum profile 

Our goal is to determine how the angular momentum 
is distributed in the halo. We wish to compute the mass 
distribution of specific angular momentum, M(< j), and 
also the spatial distribution oi j. 

We first compute the total J for each halo, and let it de- 
fine the z direction. Then we subdivide the spherical vol- 
ume of each halo into many spatial cells, each containing 
many particles, and compute the specific angular momen- 
tum in each cell, projected along the z direction. In general 
we use the symbol j to refer to this z-projected angular 
momentum. In some instances we also investigate the non- 
projected, absolute magnitude of angular momentum in 
each cell, which we refer to as When comparisons be- 
tween the two measurements are made, we explicitly add 
the subscript z to the projected value, jz- 

The cell geometry and spatial distribution were designed 
to maximize the number of particles in each cell while at 
the same time sampling as much of the j distribution as 
possible. The vector character of angular momentum nat- 
urally introduces a preferred axis, convolved with the gen- 
eral spherical symmetry of the mass distribution. The cells 
are defined using the usual spherical coordinates {r,9,(j)) 
about the halo center. Each of our cells spans the full 27r 
range in (f), and they span the range of (r/R^, sin 9) from 
(0, 0) to (1, 1). The radial shells are spaced such that each 
contains the same number of particles. The shells are then 
divided into three azimuthal cells of equal volume between 
sin0 = and 1. Positions with the same r sin 6 above and 
below the equatorial plane belong to the same cell. For 
all the M(< j) profiles studied in detail and presented in 
the figures below, we adjusted the number of radial bins 
such that each cell contains roughly Np particles, where 
Np = 500 if Mv > 5.5 x lO^^ /i-^Mg and Np = O.lM^/nip 
for Mv < 5.5 X 10^2 H-^Mq. 

We construct M(< j) profiles for each halo by ranking 
the cells by their j values, and counting the cumulative 
mass in cells with angular momentum less than j . Because 
j corresponds to a projected component, it is possible for 
a cell to have a negative value of j. The anti-alignment 
needed for a cell to have a negative j value is rare, but 
when one occurs, we remove the cell completely from the 
ranked list used to construct the M(< j) profile. About 
5% of halos have a significant amount of their total mass 

® We include all halos when we present the distribution of half-mass alignment cosines in §4.1 



(> 10%) contained in negative j cells. We do not consider 
these anti-aligned halos in our M(< j) analysis."'' 

We assign an error to the j in each cell using Eq. (7). 
Every halo analyzed has N > 1000 particles within its 
virial radius, corresponding to > 10^^ H^^Mq. Wc 
have ~ 600 halos that meet our requirements. In practice, 
we determined the shape of the functional form which fits 
M(< j) by examining the group of (~ 200) halos with 
> 6000 particles, for which the profiles are determined 
with greater accuracy. Wc extendcnl the mass range down 
to 1000 particles in order to determine whether more typ- 
ical galactic-sized halos have a similar profile. Although 
these smaller halos have larger errors on their fit parame- 
ters, they seem to obey the same M{< j) profile. 

The main limitation of this analysis is mass resolution. 
Because we need > 100 particles per cell in order to obtain 
a reasonable measure of j, our determination of M(< j) is 
onlv reUable down to ^ 10% for < 5.5 x 10^^ H-'^Mq 
and ~ (5 X 10'*mp/Mv)% for > 5.5 x lO^^ h'^MQ. 

3. MASS DISTRIBUTION OF ANGULAR MOMENTUM 

3.1. A universal angular-momentum profile 

We find that the specific angular-momentum profiles of 
halos in the simulation are well-fit by the following two 
parameter function: 



M(< j) = 



H>1. 



(8) 



30+3 

The profile has an implicit maximum specific angular mo- 
mentum jmax = io/(A* — !)• It is roughly a power-law for 
3 % 3o, and flattens out for j > jo. As can be seen by 
its relation to jmax, the quantity /i (> 1) acts as a shape 
parameter: for /x ^ 1, M{< j) is a pure power law, while 
/i — > 1 means that only half the mass falls within the 
power-law regime and the bend is pronounced. 

The pair of parameters /x and jo fully defines the angu- 
lar momentum distribution of the halo. The global spin 
parameter is related to /i and jo via 

jo bill) = V2VyR^ A' , (9) 



where 



6(m) 



-dm, 



/J, -TO 



-H ln(l-^-^)-l. 



(10) 



Once the value of A' is known, Eq. 8 is a one-parameter 
fitting function. In fitting the M{< j) profiles for each 
halo, wc first measured A' and then obtained the best fit 
M{< j) with the constraint that fj, and jo return the mea- 
sured A'. The pair (A',/i) is perhaps the most useful pair 
of parameters for characterizing the halo angular momen- 
tum because A' is the conventional global measure and jo 
can be explicitly determined from A' and /i via Eq. (9), 
while ji cannot be determined explicitly from the other 
two parameters without iteratively solving Eq. (10). Al- 
ternatively, we provide the following approximation for the 
iterative solution of Eq. (10): 



nib) ~ 1 + 



0.5 



(11) 



[exp(1.096) - 1] ' 

which is accurate to within 8% in /x — 1 over the range 
/X - 1 = 0.01- 10. 
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Fig. 4. — Mass distribution of specific angular momentum in four halos spanning a range of /i values from 1.04 to 1.9. Symbols and errors 
correspond to the ranked j measurements in cells, while the curves are the functional fits, M(< j) = MvfJ,j/{jo + j)- (a) All profiles are 
normalized to coincide at Mv, where j = j'max- The value of fi measures the relative extent of the power-law regime until it bends over. 
Shown for comparison is the distribution for a uniform sphere in solid-body rotation (dashed line), (b) All profiles are normalized to coincide 
at jo and on top of the universal profile (curve). The value of /i now correlates with the uppermost point, jmax/io> along the universal curve. 



Figure 4 shows M(< j) examples for several of our ha- 
los. In the left panel we have normalized the profiles by Mv 
and jinax in order to illustrate how different values of ^ af- 
fect the distribution. As fi approaches its minimum value 
of 1, a larger fraction of the halo mass is spinning slowly 
compared to jmax- Larger values of fi imply a more uni- 
form j distribution. The right panel shows the same four 
halos, now normalized to their best-fit jo and Mq values, 
where Mq = M{< jo) = Notice how remarkably 

the halos follow the characteristic functional shape. 

The spread in profile shapes can be described by the 
distribution of /i — 1 (= jo /jmax)- Figure 5 shows the dis- 
tribution of logiolM — 1) for all of our halos (histogram) 
along with a Gaussian distribution with the same mean 
and standard deviation: (log]^o(M ^ 1)) = —0.6, a — 0.4. 
The implied 90% range is ^ - 1 ~ 0.06 - 1.0. 

3.2. Correlations between parameters 

Although the spin parameter A' and the shape parame- 
ter /i clearly measure two different aspects of the angular- 
momentum distribution, they are not necessarily uncorre- 
lated. Figure 6 shows the joint distribution of these two 
parameters for our simulated halos. We observe a weak 
but significant linear correlation between log(/i — 1) and 
log A', with a Pearson correlation coefficient of r ~ 0.23 
(corresponding to a probability p ^ 10~® for no correla- 
tion). Thus, high-spin systems tend to have more evenly 
distributed (power-law) M{< j) profiles than low-spin sys- 
tems, but the scatter about this trend is large. 



Does the angular-momentum distribution correlate with 
the mass-density distribution? We know that the global 
spin characterized by A' does not correlate strongly with 
the halo mass (Barnes & Efstathiou 1987) and correlates 
only weakly with the concentration parameter (Fig. 2). 
Figure 7 shows n ~ 1 versus these mass parameters. We 
detect a marginal anti-correlation with mass, with r ~ 
—0.142 (p ~ 4 X 10~^ for no correlation), and an even 
less significant anti-correlation with c^, of r ~ —0.055 
{p ~ 0.26 for no correlation). 



4. SPATIAL DISTRIBUTION OF ANGULAR MOMENTUM 

In the previous section we explored the distribution by 
mass of the specific angular momentum within each halo. 
This does not tell us much about how well the angular 
momentum is aligned throughout the halo or how angular 
momentum is distributed spatially. This section is explic- 
itly devoted to these issues. In §4.1 we address the question 
of alignment, in §4.2 we investigate the cylindrical symme- 
try of the angular-momentum distribution, and in §4.3 we 
explore the angular-momentum profile in spherical shells. 



4.1. Alignment 

Here we address the question of how well the smoothed 
angular momentum is aligned throughout the halo volume. 
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Fig. 6. — The joint distribution of tiie parameters describing tlie 
angular-momentum distribution, log(/^ — 1) and log(A'), showing a 
weak correlation. The symbol size is inversely proportional to the 
relative error on — 1. The two linear regression lines are shown. 



jj.— 1 

Fig. 5. — The distribution of ^ values (histogram). Shown is a 
Gaussian distribution in log(/i — 1) (smooth curve) with the same 
mean (—0.6) and standard deviation (0.4) as the measured distri- 
bution. 

For our first alignment statistic, we divide the parti- 
cles in each halo into an inner half-mass sphere and an 
outer half-mass shell and measure the cosine of the an- 
gle between inner and outer angular-momentum vectors, 
COS01/2- The distribution of measured alignment values is 
shown in the top panel of Figure 8. 

Before we attempt to interpret these results, it is es- 
sential to keep in mind the large errors involved in deter- 
mining angular momentum direction vectors. Even for a 
halo with an intrinsically well-aligned angular momentum 
distribution (cos0i/2 — 1), measurement errors in each of 
the six directional components involved in the alignment 
cosine will tend to drive the measured value towards zero, 
and, if large enough, will lead to a false measure of mis- 
alignment. The distribution in the top panel can therefore 
serve as a lower bound for the true alignment. 

In order to obtain a conservative estimate of the distri- 
bution of true ahgnment cosines, we use Eq. (7) to assign 
errors to each component of J in both the inner and outer 
half-mass regions and perform a standard propagation of 
errors to obtain an estimate for the measurement error 
on COS01/2. We add the estimated error in each case to 
the measured alignment value in order to obtain a "cor- 
rected" distribution for cos0i/2. This is shown in the lower 
panel of Figure 8. As argued above, we believe Equa- 
tion 7 is an overestimate for the measurement error, so 
our corrected alignment cosines represent a plausible up- 
per limit for the alignment. The true intrinsic distribution 
probably lies somewhere between the corrected and uncor- 
rected distributions. If we make this assumption, then the 
two distributions allow us to conclude that between 70% 
and 90% of the halos are aligned to a greater degree than 
COS01/2 = 0.7, although there is also a tail of significantly 
misaligned halos. 

^ We also find a trend with the measured value of cos0x/2 ^i^d A'; however, this seems to be due to the correlation between J errors and the 
spin. The correlation is not apparent when we use our corrected estimate for cos6x/2- 



If the angular momentum within a halo is poorly 
aligned, one might expect a reflection of this to appear 
in its M(< jz) profile. There is a hint for such a trend in 
Fig. 8, where the distribution of cos 6*1/2 for all the halos 
is compared to the distribution for low /i halos of < 1.1. 
The latter clearly show a higher fraction of misaligned ha- 
los. Figure 9 shows the corrected cos 6*1/2 values versus 
the shape parameter ^ for each halo, and indeed there is a 
significant trend, with r ~ 0.36 [p ~ 10""'^^). ^ The mis- 
aligned halos tend to be associated with low-/i profiles — 
that is, profiles that deviate significantly from power laws, 
with a relatively larger amount of mass in the tails of their 
jz distributions. This result suggests that there are com- 
mon aspects to the origin of misalignment and that of the 
non-powerlaw nature of the M{< j) profiles. 

Since halos with misaligned angular momentum distri- 
butions may be less likely to host large disk galaxies, this 
result should be kept in mind when modeling galaxy for- 
mation within halos with low-/i angular momentum pro- 
files. 

Fig. 10 addresses the angular-momentum alignment in 
halos in another way, independent of any specific spatial 
symmetry. Here we show the M{< j) profiles for four dif- 
ferent halos, comparing the distributions of both jz and |j| 
in cells throughout the halo volume. We see that in most 
cases, the two kinds of profiles are similar; for at least half 
the mass they differ by less than a factor of two. This 
indicates that the angular momentum is reasonably well 
aligned throughout the halo. 
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Fig. 7. — Angular-momentum shape parameter versus mass 
parameters. Upper: halo mass Mv. Bottom: halo concentration pa- 
rameter Cv. Symbol sizes are inversely proportional to the relative 
errors on — 1. 
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Fig. 9. — Alignment cosine (corrected) cosSi/2 versus angular- 
momentum shape parameter ^l. The symbol size is inversely pro- 
portional to the fit error on ^L. 
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Fig. 8. — The differential distribution of the alignment statistic 
cos9x/2i measuring the alignment between the angular- momentum 
vectors within the inner and outer half mass of halos. Top: the mea- 
sured values of cos 6x/2- Bottom: after a conservative correction has 
been applied to account for the error in direction measurement. The 
true distribution should be between the distributions shown in the 
two panels. The open histograms show all halos, while the shaded 
histograms correspond to halos with /i < 1.1. The histograms are 
normalized to yield unit integrals over cos 61/2- 
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Fig. 10. — Alignment in four typical halos. Shown are the pro- 
files for the angular momentum projected in the ^-direction (filled 
points, solid line fits) and those obtained using |j| (open points, 
dashed line fits). 
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Fig. 11. — The A4(< j) shape parameters for lialos using \j\ in 
each cell, versus those obtained using jz in each cell. The symbol 
sizes are inversely proportional to the error on fi^j^ — 1. 

The slight misahgnment is expressed in a trend that the 
I J I distributions systematically resemble a power-law more 
than the jz distributions. This is demonstrated in Fig. 11, 
which shows the values of /x — 1 obtained using \j\ ver- 
sus those obtained using j^. The correlation is obvious, 
with r ~ 0.75. The value of /i — 1 is typically about a 
factor of 3 higher for the |j| profiles. The fact that this 
factor is somewhat higher for low values of /x, reflects the 
correlation seen in Fig. 9 between alignment and /U. 



4.2. Cylindrical symmetry 

Cylindrical symmetry is the natural symmetry in the 

spatial distribution of angular momentum as a result of 
tidal-torques or a sequence of mergers. In order to explore 
the degree of cylindrical symmetry, we have again divided 
the halos into cells, as described in §2.3, except that we 
no longer explicitly impose symmetry about the equato- 
rial plane and rotational symmetry in (p. Specifically, we 
define the cells as outlined in §2.3, but then, for each value 
of r sin 6, we further divide the cell particles into four re- 
gions: above and below the equatorial plane, and split by 
rotation angle = — tt and (f) = n — 27r. Because each 
cell now contains roughly one fourth of the particles, we 
increased our minimum halo mass threshold by a factor 
of four to 4 X 10^^ /i~^Mq, only for this section of our 
analysis. 

Figure 12 shows a map of the j distribution in cells in 
four of our halos. Except for the halo depicted in the 
lower-left panel, the symmetry of the angular-momentum 
distribution tends to be more cylindrical than spherical. 



with a weak dependence on at a constant distance from 
the global spin axis, r sin 9. On the other hand, the angu- 
lar momentum of the halo in the lower-left panel tends to 
be concentrated in the z = plane. 

Another way of visualizing the degree of cylindrical sym- 
metry of halos is shown in Figure 13. It shows the values of 
j in the same cells and for the same four halos displayed in 
Figure 12, now plotted as a function of the distance from 
the angular momentum axis, r/i?v sinO, and distinguished 
by different symbols according to their average distance 
from the equatorial plane, \z/Ry\. The different symbol 
types should be well-mixed if a halo is cylindrically sym- 
metric. This is the case for three out of the four halos 
shown, except in the lower-left panel, as expected from 
Figure 12. 

Although the angular momentum distributions in most 
of the halos tend to qualitatively show a certain degree 
of cylindrical symmetry, there is some indication that, for 
fixed r sin mass near the equatorial plane tends to have 
more specific angular momentum than corresponding mass 
in the poles. In order to quantify the extent to which this 
is true, we compare the quantity jz=o/jpo\c, defined as 
the ratio of specific angular momentum in cells near the 
equatorial plane {\z/Ry\ = — 0.35) to that in cells near 
the poles {\z/Rv\ = 0.65 — 1). For each halo, we deter- 
mine jz=o/jpoie by averaging over all values of rsin^ for 
which the two z ranges overlap. For example, in Figure 13, 
Jz=o/ipoic can be estimated for each halo by comparing the 
average value of j for the solid squares {\z/Rv\ ~ 0.17) to 
that for the open circles {\z/Ry\ ~ 0.83) at fixed values of 
r sin 9. Figure 14 shows the histogram of this ratio for all 
of our halos. The mean is (log]^oOz=o/jpoio)) — 0.13, and 
the standard deviation is 0.18. 

In order to test whether our specific choice of cell ge- 
ometry would bias this measure, we also generated 500 
realizations of halos with perfect cylindrical symmetry 
Oz=o/jpoie = !)• They each had NFW density profiles and 
M{< j) profiles in the form of Equation 8. The resulting 
jz=o/ipoie histogram for the artificial catalog is shown by 
the shaded historgram in Figure 14. The sharp peak near 
a ratio of 1 indicates that the bias due to our cell geometry 
is much smaller than the detected deviation from cylindri- 
cal symmetry. Recall however that the average deviation 
from cylindrical symmetry is small, about 35%. 




Fig. 12. — Map of the angular momentum distribution in four representative halos. The shading code is shown, with the j'max cell in each 
halo shaded dark and the minimum j in each halo shaded light. The arrow indicates the direction of total J in each halo. A higher quality 
color .eps figure is provided seperately on astro-ph. 



A UNIVERSAL ANGULAR MOMENTUM PROFILE FOR GALACTIC HALOS 



11 




: j(r)«M 



0.1 



0.01 



0.001 




(r/R^,i,,)sine 

Fig. 13. — Cylindrical symmetry? Each panel shows the spe- 
cific angular momentum in spatial cells for the same halos shown in 
Fig. 12. The value of j/jmax for each cell is plotted as a function 
of the distance from the angular momentum axis, rsinS. The three 
symbol types represent the cell's average distance z/R^i^ from the 
(z = 0) equatorial plane, as indicated in the figure. 
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Fig. 15. — Specific angular-momentum profiles in spherical shells 
as a function of the mass encompassed by that shell. Shown are 
profiles for three halos from our simulation (symbols), and the cor- 
responding power-law fits, normalized at the outermost shell. 
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Fig. 14. — Deviation from cylindrical symmetry. Histogram of 
jz=o/jpole for all our halos (solid line), defined as the ratio of the 
specific angular momentum about the 2 = equatorial plane and 
that near the poles, averaged at fixed distances from the angular 
momentum axis. The bias due to the specific choice of cells is cal- 
ibrated by the distribution of this statistic for an artificial set of 
purely cylindrically symmetric halos (shaded). The histograms are 
normalized so that they have unit integrals over logiQ(jrz=o/ipole)- 
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Fig. 16. — The slope s of the angular-momentum profile in spher- 
ical shells, from the fit j{r) oc M", versus halo mass. The symbol 
size is inversely proportional to the fit error on s. The linear re- 
gression lines are shown (thin lines). Shown as a thick line is the 
prediction from linear tidal-torque theory (§5). 
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Fig. 17. — The slope s of the angular-momentum profile in spher- 
ical shells versus the shape parameter ^ of the mass distribution of 
angular momentum, A/{< j). The symbol size is inversely propor- 
tional to the fit error on /i — 1. Shown are the two regression lines. 



4.3. Profile in spherical shells 

Although, as demonstrated in §4.2, the specific angular 
momentum distribution in halos is closer to being cylin- 
drically symmetric than spherically symmetric, an explo- 
ration of how j, averaged over shells, behaves as a function 
of radius is of particular interest because this is the profile 
that is most directly addressed by theory (see below). 

Barnes & Efstathiou (1987) found in their early work an 
indication for a rough j(r) cx r behavior. For our sample 
of well-resolved halos, we indeed find that a power-law: 
j(r) oc r", with a roughly distributed over the halos like 
a Gaussian: a = 1.1 ± 0.3 is a good approximation. 

In order to connect more closely to the mass distribution 
profiles discussed in previous sections, and also to allow a 
more direct connection with theory, we focus on charac- 
terizing the spherical profile as a function of the mass con- 
tained within the radius r, M — M{< r). We find that the 
power-law approximation describes the profiles reasonably 
well, 

j(M) cx iVr , (12) 

with 8 roughly distributed over the halos like a Gaussian: 
s = 1.3 ±0.3. 

Figure 15 shows three examples of the spherical j{M) 
profiles along with their best-fit power-laws. Notice that 
although the power-laws serve as a good general character- 
ization of the j behavior with M, the profiles do not always 



follow a power-law form in detail. In many cases, the pro- 
files show an upwards bend in the outer shells, reminiscent 
of the characteristic bend in the low-// mass distribution, 
M{< j). In some cases, j is not always monotonic with 
M. So, although the s values provide a useful character- 
istic of the spherical distribution of j , some information is 
lost under this approximation. 

Shown in Fig. 16 is the distribution of s values, ver- 
sus halo mass. Despite the large scatter, there is an 
anti-correlation between slope and mass, with r ~ —0.49, 
{p ^ 10~^^). As we will discuss in §5, a simple calcula- 
tion based on linear tidal-torque theory provides a hint at 
understanding the typical slope value s > 1, and predicts 
a mild anti-correlation between s and the halo mass, in 
agreement with the simulations. 

The spherical profile j{M) and the mass distribution 
M{< j) are expected to be related. If the mass distribu- 
tion is spherically symmetric and the angular-momentum 
distribution is cylindrically symmetric, this relation can be 
spelled out exphcitly. In this case the inverse of M(< j) is 
the spatial profile jcyi{M) in growing cylinders about the 
total spin axis. A power law with a bend for M(< j) trans- 
lates to a corresponding power law with a weaker bend in 
the opposite direction for j{M). We thus expect a corre- 
lation between the shape parameter n of M{< j) and the 
slope s of the spherical profile j{M). A high s and a low 
/i both correspond to halos with a high fraction of mass 
spinning slowly relative to jmax- Figure 17 shows s and /i 
for our halos. An anti-correlation is evident, as expected, 
with r ~ —0.36 [p ~ 10"^'^). The scatter refiects devia- 
tions from the global spatial symmetries. It is likely, for 
example, that partial misalignment and pockets of low- 
angular momentum material at large radii contribute to 
the scatter. By comparing Figures 16 and 7 it is evident 
that the scatter in the s-/i relation is large enough to wash 
away any mass dependence of fx that would be implied by 
the weak mass dependence of s. 

5. ON THE ORIGIN OF THE PROFILE 

In this section we explore the possible origin of the an- 
gular momentum profile based on two different pictures 
for how specific angular momentum is acquired in halos. 
In the first picture, angular momentum is built up shell by 
shell, and is modeled using linear theory with no further 
angular momentum transfer within the halo. In the sec- 
ond picture, the angular momentum profile arises solely 
from angular momentum transfer from one or more satel- 
lite halo merger events. Since halo mass accretion histo- 
ries typically refiect some combination of relatively quies- 
cent mass accretion as well as more pronounced mergers 
(where j transfer is likely) , the j profiles of individual halos 
may refiect a complicated combination of the two processes 
explored here. Indeed the two processes are intrinsically 
linked, since, in the hierarchical framework, all mass ac- 
cretion can be treated as the accretion of smaller objects. 
As shown below, the profiles calculated from linear theory 
as well as those calculated using j transfer from mergers 
produce a range of /i values, in qualitative agreement with 
what is observed in the simulations. 

5.1. TTT and EPS 



A UNIVERSAL ANGULAR MOMENTUM PROFILE FOR GALACTIC HALOS 



13 



Linear tidal-torque theory (TTT) provides a hint for the 
origin of the detected power-law-hke j profile, as follows. 
TTT (Doroshkevich 1970; White 1984) implies that the 
angular momentum gained by a halo at time t before its 
turn-around is 



Equation 14 provides the required value of j associated 
with each mass shell given its "accretion time", tc{M). 
Considering, for example. My and Mv/2, we obtain 

C7(Mv) 



s(Mv) -1.44 In 2^/3 £)[ti/2(^^)] 



J,{t)=aitfD{t)eijkTjiIik , 



(13) 



-3/2. 



where the time growth is from some fiducial initial time, 
a{t) is the expansion factor at t, D{t) is the linear growth 
factor, Ilk is the inertia tensor of the proto-halo at the 
initial time and Tji is the tidal tensor at the halo center 
at the initial time, smoothed on the halo scale. This is 
based on assuming the Zel'dovich approximation for the 
velocities inside the proto-halo, and a 2"'^-order Taylor ex- 
pansion of the potential. We show elsewhere (Porciani, 
Dekel & Hoffman 2000) that the standard scaling relation 
of TTT should be slightly modified; it should read 



(14) 



a{M,/2)J ' ^^^^ 
where ti/2 is the time at which the halo of mass 
first accreted half of its mass, ti/2{My) = tc{My/2). We 
adopted for all halos, by the requirement of virialization, 
tc(M^.) = to (the current age of the universe), and set 
D{to) — 1. The most probable value of ti/2{My) can be 
estimated via EPS (Lacey & Cole 1993, Eq. 2.26). The 
predicted s{M^) relation is added to Fig. 16 as the thick 
solid line. The agreement with the mean trend shown 
by the simulated halos is remarkable, especially since our 
model has no free parameters. 

We can further elaborate the model by following the 
more detailed mass accretion history of halos using EPS 
merger trees. Specifically, we have used the method of 
Somerville & Kolatt (1999) to model mass growth and halo 
accretion histories for a random ensemble of dark matter 
halos formed in the ACDM model and with masses in the 
mass range studied in the simulation. At each epoch, ti, 
and time interval, AU, the model draws a mass AMi ac- 
creted by a halo of mass Mi{ti) using the EPS probability 
distribution. If we assume that this mass is accreted with 
specific angular momentum given by Eq. (14) and that 
there is no angular momentum loss, we can integrate the 
specific halo mass accretion history to get the M(< j) 
distribution for the final halo at 2; = 0. Figure 18 shows 
M(< j) profiles calculated using this model along with 
the fits of the form M(< j) = M^nj{j + jo)"\ Eq. (8), 
for three representative mass accretion histories from the 
ensemble of realizations. One can see that the model re- 
produces the general shape of the M(< j) profile with 
values of that are similar to those of simulated halos. 

For a given halo mass, the differences in merger histo- 
ries lead to scatter in the resulting values of n, with the 
range of — 1 ~ 0.1 — 2.5. The absence of small fj. (< 1.1) 
may result from our simplified model assumption of perfect 
alignment. Indeed, the range of these modeled /i values is 
similar to the range obtained for our simulated halos when 
\j\ profiles are considered instead of jz (see Fig. 11). 

We can also reproduce both the range of values and the 
trend with mass of the slope s if we force a power-law fit 
to the model j{M) profiles. The halos with low (high) /i 
in this model are those that accrete large fractions of their 
mass at early (late) epochs. 

Interestingly, the value of /x for a given halo in this model 
strongly depends on the halo's formation redshift, Zf, de- 
fined as the redshift where the mass of the most massive 
progenitor was half of the final mass for the first time 
[n — 1 oc {1 + Zf)~'^]. In light of this expectation, it is 
somewhat surprising that we find no noticeable trend with 
redshift for ji parameters in halos of fixed mass (see Ap- 
pendix) but we are currently investigating whether a trend 
is evident using individual halo merger histories (Wechsler 
et al. 2001). If there is a trend with formation time, it 
ln(M2/Mi) ' ^""^ may have interesting implications for galaxy formation. 

The j dependence on the time of collapse is oc tc rather than oc t^^^ , because the tidal part of the deformation tensor, the source of angular 
momentum, is actually independent of the initial density of the perturbation, which determines the collapse time. The shear depends instead 
on the rms fluctuation <t, which involves an implicit mass dependence. 

® The time of accretion is expected to be proportional to the turnaround time; see e.g. Sugerman et al. (2000). 



where j is the specific angular momentum of the halo, tc 
is the turn-around time of the halo, and a{M) is the rms 
density fluctuation on scale M at the initial time. If the 
perturbation turned around while the cosmology was still 
Einstein-de Sitter, then j octc- 

To understand the general power-law behavior of the 
angular- momentum profile we apply Eq. (14). We can 
assume that mass is accreted in shells and that the turn- 
around_time of each shell is determined by 5{M)D{tc) ^ 1, 
where 5{M) is the mean density inside M . For a Gaussian 
field, the typical density fluctuation profile about a ran- 
dom point scales like (Dekel 1981) 5{r) oc £,{r)^ where ^(r) 
is the linear two-point correlation function. This is accu- 
rate to a few percent also around a high peak (Bardeen et 
al. 1996, Fig. 8). Thus, for a power- law power spectrum 
Pfe a fc" and a flat universe, we obtain the following spin 
profile within each halo: 

i(M) oc m2/3+(3+")/3 . (15) 

This implies 1 < s < 4/3 for — 2 < n < —1, the range ap- 
propriate for 10^^ — lO^^/i~^M0 halos in a CDM spectrum, 
in pleasant agreement with our finding for the simulated 
halos (Fig. 16). 

The tidal torque theory can be combined with the Ex- 
tended Press-Schechter formalism (EPS, Bond et al. 1991; 
Lacey & Cole 1993) to obtain a more detailed model of j- 
distribution, as well as quantitative predictions for values 
of s and ji and their mass dependence. 

Let us assume for simplicity that halo mass is acquired 
via accretion of material with specific angular momentum 
given by Equation 14 and that the direction of angular 
momentum of each mass shell is perfectly aligned. 

We can readily estimate the slope s of the j{M) oc 
distribution if we assume that the turn-around time for 
each mass shell is approximated by the time at which the 
mass was first accreted onto the halo^ by comparing the 
expected j value at any two shells encompassing masses 
Ml and M2, 

_ln[j{M2)/j{M,)] ^^g^ 
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Fig. 18. — Angular-momentum profiles calculated using linear 
theory and extended Press-Schechter merger trees. The jagged lines 
are the derived profiles and the smooth lines are the best fit profiles 
using Eq. (8). 



5.2. Profile due to minor mergers 

A simple toy calculation for how angular momentum 
is being deposited in a large halo during a minor merger 
seems to provide another hint for the origin of the charac- 
teristic j profile. This is studied in more detail using toy 
models and full A^-body simulations by Dekel & Burkert 
(2000). 

Consider a fixed halo and an incoming satellite halo of 
mass profiles M{r) and m{i) respectively. Assume that 
as it moves in, the satellite is losing mass outside a tidal 
radius it , which we approximate at r to be determined by 



m(4) _ itn{r) 



K(r) = 



2M{r) - r 



dM 



(18) 



For host halos and satellite halos with power-law mass pro- 
files (M oc r", a < 2), the above approximation implies 
that the total mass of a satellite halo located at a radius 
r of its host is m(r) = m[£t{r)] oc M(r). 

Suppose that the satellite is spiraling in due to dynami- 
cal friction roughly along circular orbits, and that the mass 
and J are deposited locally. Then the spherical j{r) profile 
is obtained by averaging over shells. 



Anr'^ p{r) j (r) ~ m{r) 



d[rVc(r)] dm(r) 



dr 



dr 



rVcir) . (19) 



The two terms on the right hand side reflect angular- 
momentum transfer due to the slowdown of the satellite by 
dynamical friction and the direct tidal stripping of mass 
from the satellite halo, respectively. For an isothermal 
host and satellite, we obtain 



..0.1 



m(r) = consL. 
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Fig. 19. — Angular-momentum profile in spherical shells due 
to a minor merger of a satellite into a fixed NFW halo, assuming 
circular orbits and three difi'erent recipes for tidal mass stripping. 
The marked heavy solid lines are for mass-loss recipe m(r) oc M{r) 
and for the mass loss in the case of an isothermal satellite profile 
respectively. The thin solid line is for the limiting case of no mass 
loss. 

Fig. 19 shows the expected j{M) profile based on this 
toy model for an NFW halo swallowing a satellite, using 
three alternative recipes for mass loss by the satellite: (a) 
m(r) oc M(r), (b) a satellite with an isothermal profile, 
and (c) a satellite of fixed mass, for comparison. The gen- 
eral shape of j(M(< r)) in the outer decade of mass, when 
realistic mass loss (a or b) is considered, is reminiscent of 
the profiles detected in the cosmological simulation (see 
Fig. 15). ^ ^ 

Full A"-body simulations of mergers, spanning a range 
of halo and collision parameters, confirm the robust pro- 
duction of such characteristic profiles (Dekel & Burkert 
2000). The resulting j{M) from a sequence of minor merg- 
ers is expected to be a sum of similar contributions, each 
projected onto the direction of the total net angular mo- 
mentum (perhaps determined by the most major merger 
involving a large fraction of the final mass). 



6. DISK STRUCTURE 

We now explore certain possible implications of our re- 
sults on the angular-momentum distribution in the halos 
for the formation of galactic disks. We follow the lim- 
iting case assumption that the specific angular momen- 
tum is conserved during gas infall into a centrifugally sup- 
ported disk, and compare the resultant gas density profile 
to an exponential disk. Also, since the angular momentum 
structure of halos is governed by at least two parameters 
— say ^, in addition to A' — we specifically explore how 
variations in the shape parameter may affect the disk char- 
acteristics. 
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Fig. 20. — The disk surface-density profile (heavy solid line) 
implied by our M{< j) distribution with the indicated parameters 
under the assumption that specific angular momentum is conserved 
during adiabatic baryonic infall. Shown for comparison are an expo- 
nential disk and the disk profile resulting under similar assumptions 
from a uniform, solid-body rotating sphere. 

Suppose that each dark-matter halo has a specific an- 
gular momentum profile M{< j) given by Eq. (8). Given 
fi, M{< j) is determined by the spin parameter via 
Jo = \/2\'VvRvb{fj,)^^ . Assume that initially the specific 
angular momentum distribution of the gas mirrors that of 
the halo, Mgas(< j) — fAU{< j), where / is the mass 
fraction of the halo that ends up as cool baryons in the 
disk, Mdisk = /Mv. 

After infall into a disk while preserving j, the specific 
angular momentum of a gas element that ends up in a cir- 
cular orbit of radius r is j{r) ~ -\/ GM {r)r, where M(r) is 
the total mass of DM and baryons within r. All the gas 
with specific angular momentum less than j{r) will wind 
up interior to Using Eq. (8), we obtain for the disk mass 
profile 

md(r) ~ //iAfv . ^^''l , , Jir) < J„.ax. (21) 

In order to estimate the implied disk structure, and to 
gain a qualitative understanding of how the value of /i 
may affect the disk size and density profile, we assume 
first that the dark halo does not react to the infall of gas 
and approximate the total mass distribution of the system 
to be that of an isothermal sphere, M{r) cx r. In this case 
j(r) ~ rV{r) = rV^. The final mass distribution of the 
disk is then 

r 



r <r„ 



Here, rd = V^X' Rvb{fi)~^ and r„ 
implied surface density profile is 



(22) 



rd/{^l - 1). The 
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In our simplified model, the surface density vanishes be- 
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Fig. 21. — The resultant disk surface-density profiles, as in 
Fig. 20, for three different values of fi. The abrupt edges of the 
disks correspond to j = jmax in each case. Such an edge is clearly 
unphysical; they arise from our simplified functional characteriza- 
tion of the A/(< j) profiles of halos. True profiles are likely to have 
smoother edges, the study of which requires simulations with higher 
mass-resolution. 



Several implications are apparent from the derived disk 
surface-density profile of Eq. (23). First, the distribution 
is qualitatively similar to an exponential disk — it falls 
rapidly at large r, and flattens towards the origin. How- 
ever, in detail, there are significant differences. At large 
radii, r ^ r^, the fall off Tjd{r) oc is slower than expo- 
nential. Near the origin, r <C r^, the profile is steeper than 
exponential, Sd(r) cx versus cx const. As expected, the 
scale radius scales linearly with A'. However, the scale 
radius is also a strongly decreasing function of /z. For a 
fixed spin parameter A', as fi varies over its ~ 90% range 
about the mean (/i — 1 = 0.06 — 1.0), rd varies by about 
a factor of 5. This is similar to the variation of rd as a 
function of A' for fixed /i (a factor of 5 over the 90% range 
of A'). 

Figure 20 shows the disk surface-density profile derived 
using a more accurate numerical calculation. Here we as- 
sume that the dark halo initially followed an NEW density 
profile with Cv = 14. The M(< j) profile is characterized 
by A' = 0.04 and fi = 1.25, and the disk mass fraction is 
assumed to be / = 0.03. The disk profile is obtained by 
self-consistently solving for the disk circular velocity at r 
and taking into account the adiabatic contraction of the 
halo in response to the baryonic infall (Blumenthal et al. 
1986). Shown in comparison is the surface density pro- 
file obtained under the same assumptions regarding the 
halo structure and j conservation except that the given 
M(< j) distribution is replaced with the assumption that 
the final disk profile is exponential. A similar result is ob- 
tained by assuming that the initial M(< j) is that of a 
uniform-density sphere in solid body rotation. 

Figure 21 shows how the derived disk profile varies as 
a function of /i. As expected, the central concentration of 
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the disk is a decreasing function of fx. 

As anticipated, compared to the exponential disk (and 
solid body rotator) the derived disk profiles have core and 
tail excesses. Only halos with the highest /i values (/i > 2, 
corresponding to ^ 5% of the halos) produce surface den- 
sity profiles that begin to resemble an exponential disk. 

In comparison with the observed light profiles of disks, 
the derived gas surface density profiles deviate in their gen- 
eral shape and in particular in their central concentration. 
The general deviation of shape is less worrisome because 
the final stellar profile may not necessarily mirror the gas 
profile. As Lin & Pringle (1987) and Olivier et al. (1991) 
argued, viscous transport processes can act to produce ex- 
ponential stellar disks under rather general assumptions 
about the nature of the initial gas disk. However, since 
typical viscous processes would tend to amplify the central 
density even further, the indicated central density excess 
is of some concern. 

This result is reminiscent of the more severe problem of 
excessively small disk sizes seen in hydrodynamic cosmo- 
logical simulations (e.g., Navarro & Steinmetz 1997, 2000). 
That we find a disk central-density excess even under the 
assumption that angular momentum is conserved during 
the collapse indicates that the solution to the problem 
may not be trivial, because any simple angular-momentum 
transport mechanism would drive the angular-momentum 
from the gas to the halo and from the inside out and could 
only worsen the discrepancy. The solution may require 
that the baryons somehow obtain an angular momentum 
distribution that is biased relative to that of the dark mat- 
ter towards high-j. 

Alternatively, if the resultant profiles are interpreted to 
represent both a bulge and a disk component, with the 
slowly rotating core associated with the bulge, then the 
central concentration problem may be resolved. Testing 
this hypothesis will require more detailed modeling. In 
particular, it will be necessary to match the observed dis- 
tribution of galaxy bulge-to-disk ratios as well as disk scale 
lengths. 

Another way out may be associated with the fact that 
some of the approximations adopted here may break down 
when the j distribution is significantly misaligned through- 
out the halo, which indeed preferentially happens for low-fj. 
halos (§4.1). 

7. CONCLUSIONS 

We have studied the angular momentum structure 
within dark-matter halos in the mass range 10^^ — 
10^^ /i~^Mq, using a sample of ~ 600 halos in an N-body 
simulation of the ACDM cosmology. We examined in de- 
tail the more massive of these halos, for which the angular 
momentum is measured more accurately, and then verified 
that the more typical galaxy-sized halos have a similar an- 
gular momentum distribution. 

Our primary result is that the mass distribution of spe- 
cific angular momentum in halos obeys a roughly univer- 
sal form, which is well fit by the two-parameter function 
M(< j) = iJ.Mvj/{jo + i). For a fixed value of the spin 
parameter A', the parameter n determines the shape of the 
distribution, with high-^ corresponding to a pure power- 
law and a smaller contrast between the tails of the distri- 
bution. The distribution of /x — 1 is roughly log normal. 



with the 90% range spanning ji— 1 = 0.06—1.0. Compared 
to an idealized uniform sphere in solid-body rotation, the 
simulated halos tend to have more of their mass in the tails 
of the distribution, especially at small j. The shape pa- 
rameter is only weakly correlated with the spin parameter, 
fx—loc {X'Y with (5 = 2 ± 1.5 and a correlation coefficient 
r ~ 0.23. 

Most halos have well-aligned angular momentum vec- 
tors throughout their volume. Between 70% and 90% have 
alignment cosines greater than 0.7 between the inner and 
outer half masses. However, at least 10% of the halos 
have significant misalignments, which may be interesting 
for several reasons. First, they tend also to be halos with 
low /i, indicating that the direction of J in accreted ma- 
terial, or the halo merger history, plays a role in deter- 
mining the shape of the M(< j) profile. In addition, it 
is unlikely that the standard simple picture of disk for- 
mation could be valid within a severely misaligned halo, 
so perhaps these objects typically host spheroidal stellar 
components. Marginally misaligned halos may play a role 
in the formation of galactic warps (e.g., Dekel & Schlos- 
man 1983). 

The spatial distribution of halo angular momentum 
tends to be more cylindrically symmetric than spherically 
symmetric. At a fixed distance from the rotation axis, 
mass near the equatorial plane typically has only about 
~ 35% more specific angular momentum than correspond- 
ing mass near the poles. 

The mean j is spherical shells encompassing mass M is 
well-fit by a power law, j oc M*. The power s is distributed 
like a Gaussian with a mean of s = 1.3 and standard de- 
viation a = 0.3. 

We pointed out two possible explanations for the origin 
of the universal angular-momentum profile. The first is 
based on applying a corrected version of the linear tidal- 
torque theory (Porciani & Dekel 2000) to extended Press- 
Schechter mass accretion histories. The second is based on 
the nonlinear process of angular momentum transfer from 
satellite orbits in halo mergers (see also Dekel & Burkert 
2000). Since halo mass growth arises as a combination of 
merger events and relatively quiescent mass accretion, the 
origin of the profile in individual halos may reflect some 
combination of the two processes. Each of these processes 
seems to produce j profiles that are similar in form to 
those observed in simulated halos, and thus may provide 
a starting point for a deeper understanding of the origin 
of halo angular-momentum structure. 

Finally, we have started to explore the implications of 
our universal M{< j) distribution in the context of the 
archetypical model of galactic disk formation, namely, adi- 
abatic baryonic infall to a rotationally supported disk 
while conserving angular momentum in every mass ele- 
ment. We find that the implied surface-density profiles, 
which vary as a function of both fi and A', deviate signifi- 
cantly from an exponential disk for all but the largest val- 
ues of fi. The resultant surface profiles are more extended 
than exponential at large r, and are overly concentrated at 
small r. Since the observed light profiles of disks are closer 
to exponential, our result indicates that the disk formation 
process cannot be fully understood within this simplified 
picture. 

We mentioned that the general deviation from an expo- 
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nential surface density at moderate radii might be over- 
come in this picture if a viscous transport mechanism is 
included (Lin & Pringle 1987). As shown by Ohvier et 
ah (1991), such a mechanism tends to produce exponen- 
tial stellar distributions regardless of the initial gas profile. 
However, most of the expected j transport typically acts 
to increase the central mass density, so the deviation at 
small radii may be a significant problem. Even though 
we have assumed that the baryons experience no angular- 
monicintuni loss, wc ended up with an "angular momentum 
problem" similar to that detected in hydrodynamic sim- 
ulations (e.g., Navarro & Steinmetz 1997, 2000, where it 
is associated with j transport from the gas to the dark 
matter and from inside out). This problem is most se- 
vere for halos with small /i, so the tendency for these ha- 
los to have misaligned angular- momentum distributions, 
which makes them less-likely hosts of large disk galaxies, 
may partly case the problem in the limiting case where 
angular-momentum is conserved. However, even if we fo- 
cus on halos that tend to be well aligned, with > 1.1, the 
central densities remain higher than those in exponential 
disks. 

A solution to this problem may be to associate the de- 
rived central mass concentrations with bulges. It is indeed 
possible that the central regions of disks with profiles as in 
Fig. 21 are unstable to self gravity. More detailed modeling 
is needed for testing the viability of this solution. In partic- 
ular, one may be worried about the model over-predicting 
bulge to disk ratios and in particular under-predicting the 
number of bulge-less galaxies. On the other hand, such an 
interpretation may shed new light on the very origin of the 
Hubble sequence, since both /j, and A' may play a role in 
determining disk scale lengths, surface brightnesses, and 
bulge fractions. 
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APPENDIX 

Since our error estimates are only approximate, it is important to carry an independent test of the validity of our 
M(< j) profiles in view of the finite mass resolution of the simulation. The test utilizes the output from another ACDM 
simulation with identical cosmological parameters and the same number of particles as used in our main simulation, but 
now in a box of side 30 h~^Tslpc, namely half the original size. The mass resolution is thus 8 times higher, so a statistical 
comparison with halos analyzed from our lower resolution simulation at fixed mass will provide a useful test for the effects 
of mass resolution. The smaller box simulation also allows us to extend our analysis down to much smaller galaxy-mass 
halos (Mv ~ 10^^ h~^MQ). Unfortunately, the higher resolution simulation was stopped at 2; = 1.7, so we cannot compare 
directly to the z = results presented in the main paper. However, a comparison performed at high z can perfectly serve 
our purpose. This also allows us the opportunity to check for any evolution in the M(< j) profiles with rcdshift. 

Figure 22 shows the best-fit /x parameters as a function of halo mass for halos from our main simulation at 2: = 3. The 
distribution of values is very similar to that seen at 2; = (Fig. 7), so there is no indication of evolution. Figure 22 also 
shows /i versus My for halos from the higher resolution simulation at the same time. Although the overlapping region in 
mass is rather small, there is no indication of any offset between the two simulations. The halos have a similar distribution 
of shape parameters over about three orders of magnitude in mass, ~ 10^^ — 10^^ H^^Mq. 

We conclude that mass resolution does not seem to limit our analysis, and also that the M{< j) profiles of halos show 
very little variation as a function of redshift and mass. 
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Fig. 22. — The effect of mass resolution: /j, versus Mv (at z = 3) for halos in our main simulation (solid symbols) and in a simulation with 
8 times higher mass resolution (open symbols). The symbol size is inversely proportional to the relative errors on /i — 1. 
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